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Abstract 

We analyze the evolution of energy density fluctuations in cosmological 
scenarios with a mixture of cold dark matter and quintessence, in which 
the quintessence field is modeled by a constant equation of state. We ob- 
tain analytic expressions for the time evolution of the quintessence pertur- 
bations in models with light fields. The fluctuations behave analogously 
to a driven harmonic oscillator, where the driving term arises from the 
inhomogeneities in the surrounding cosmological fluid. We demonstrate 
that the homogeneous solution, determined by the initial conditions, is 
completely sub-dominant to the inhomogeneous solution for physically 
realistic scenarios. Thus we show that the cosmic microwave background 
anisotropy predicted for such models is highly insensitive to the initial 
conditions in the quintessence field. 

1 Introduction 

Growing observational evidence suggests that the total matter density of the 
universe is significantly less than the critical density. Yet, measurements of the 
cosmic microwave background (CMB) indicate that the universe is flat and the 
total energy density is precisely equal to the critical density 0. Candidates 
for the missing energy are the cosmological constant (A) [pj, and quintessence 
(Q), a time- varying, spatially inhomogeneous component with negative pressure 
H f|. Examples of quintessence include slowly evolving fields or topological 
defects such as a network of light and tangled cosmic strings §-JlJ|. Cold 
Dark Matter (CDM) models with quintessence (QCDM) which fit the data from 
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observations of high red shift supernovas, gravitational lensing, CMB anisotropy, 
and structure formation have been found. 

For the purposes of this paper, we model quintessence as a scalar field 
quintessence rolling down a potential V(Q) with an equation of state w = 
Pq/pq, where pq is the pressure and pq is the energy density of the field. We 
consider here a large class of models in which the field and fluctuation evolution 
may be represented by a constant equation of state with value between and 
-1, and in which the sound speed in quintessence fluctuations, c^q, approaches 
unity at scales much smaller than the horizon. Most models of quintessence that 
have appeared in the literature satisfy these conditions. This class includes not 
only "tracker models", which have dynamical attractor behavior, but also more 
general potentials H-jy|. We show that, in such models, the CMB anisotropy 
is insensitive to initial conditions of the quintessence field. By insensitive we 
mean that the fractional change in the CMB anisotropy power spectrum due to 
a change in initial conditions cannot be observationally resolved. 

These conclusions were originally mentioned in our first paper on quintessence 
models j|, and also by others || |[ who noted that the amplitude of the 
energy density perturbations in the background fluid are largely independent of 
the initial conditions in the scalar field, provided the initial energy contrast is 
less than unity. Our present work is consistent with these earlier conclusions. 
However, we go further in this paper by examining in detail the behavior of the 
scalar field perturbations, and the causes for the insensitivity, as described next. 

In the next section we derive the equations of motion for the evolution of 
the Q-field and its fluctuations. We introduce a formalism to numerically study 
the quintessence fluctuations in terms of the evolution of the equation of state 
as a function of cosmological scale-factor. In section || we show that models 
with constant equation of state represent the behavior of a large class of models 
with light fields and monotonically evolving equations of state. In section |] 
we analytically solve the fluctuation equation for constant equation of state at 
large wavelengths. We use this solution to describe key features of the fluctua- 
tion evolution obtained from the numerical integration. Finally, in section || we 
report the effect of changing initial conditions on quintessence fluctuation evo- 
lution and CMB anisotropy. We show that the anisotropy is highly insensitive 
to changes in the initial conditions. 

2 Field and fluctuation equations 

We consider a matter-quintessence Lagrangian of the form 

L = L B + L Q , (1) 

where the B refers to all the background species of particles and fields, including 
baryons, photons, cold dark matter, and neutrinos. The background cosmology 
is described by the Friedman-Robertson- Walker metric with positive signature. 
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We model the Q-field as a classical, self-interacting, scalar field, minimally cou- 
pled to other constituents of the universe through gravity, 

Lq = \q^Q»-V{Q). (2) 

We only consider models with canonical kinetic energy terms, since c 2 q ap- 
proaches unity at scales much smaller than the horizon Q in these models. 

The time evolution of the quintessence field is determined by the equation 
of motion 

Q" + 2-Q' + a 2 ^ = 0, (3) 

where a is the scale-factor, normalized to unity today, and where the prime (') 
represents d/dr, the derivative with respect to conformal time. By specifying 
the functional form of the potential, V(Q), along with initial conditions Q, Q' 
at a time Ti n u , the subsequent evolution is determined for all times r > Ti n u ■ 

At the scales of interest for the cosmic microwave background (> IMpc), 
fluctuation amplitudes of the Q-field and the metric are usually small compared 
to the fields themselves, and thus a linearization of the Q-field and Einstein 
equations in the perturbation suffices to describe the fluctuation dynamics. Our 
analysis, carried out in the synchronous gauge, uses the conventions and equa- 
tions from Ma and Bertschinger Jl7[ ] . 

The synchronous gauge is defined by the condition that the time-time and 
time-space part of the metric are not perturbed. The perturbed metric is given 
as: 

ds 2 = a 2 (r){-dT 2 + (Sij + h^dx^x 3 } . (4) 

The metric perturbation hij can be decomposed into a trace part h = ha and a 
traceless part r\. To enable us to study fluctuations as a function of wavelength, 
we will work in Fourier space (k) with perturbations h(k,r) and r](k, r). 

The time evolution of the perturbed metric is obtained by linearizing the 
Einstein Equations in Fourier space (see ]l7| ). The cosmological perturbation 
equations for the dynamical Q-component are obtained by expanding the scalar 
field equations about the homogeneous background. In the synchronous gauge, 
we can write down the equation for small fluctuations SQ in Fourier space: 

5Q" + 2^SQ' + (a 2 ^ + k 2 )8Q = ~h{k, t)'Q'. (5) 

The fluctuation amplitude evolves in time like a scalar field. 

To compute the quintessence fluctuation evolution and CMB anisotropy 
power spectrum, we use the fluctuation and Einstein equations to evolve the 
moments of the photon distribution in a Boltzmann code. We modified two 
separate Boltzmann codes, CMBFAST (||) and LINGER (@), by adding in 
scalar field evolution, and repeated all computations in the conformal gauge. 
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The results from these multiple approaches were all identical to better than 1 
part in 10 8 @. 

We can parameterize quintessence models in terms of the evolution of the 
equation of state as a function of scale-factor w(a). The formulation is useful 
in studying the time dependence of quintessence fluctuations and in unearthing 
the similarities in time evolution of models with different potentials. The energy 
density and pressure of the Q-component in this parameterization are given by 
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where H is the present-day value of the Hubble constant, and £Iq is the present- 
day Q-field energy density as a fraction of the critical energy density. 

To obtain the equations of motion in terms of the equation of state, we 
express the first and second derivative of the potential in terms of w: 
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If we define 5tp = 6Q/y/l + wM and ip' = Q' / y/T+wfr), we can use Eqs. 
H and H to convert the field (Eq. and fluctuation (Eq. [j]) equations for the 
Q-field into a form in which the potential is implicit: 



<•" ((1 + 3^)^)^=0, 
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A practical consequence of this change of variables is that w" drops out of the 
evolution equation, and we only need to specify w and w' as functions of time. 
Thus, given a well sampled table of the equation of state history w(a), we can 
numerically integrate this equation in a Boltzmann code to obtain the evolution 
of the quintessence field and its fluctuations. 
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3 Why constant w approximates well a large 
range of potentials 



For quintessence models described by a scalar field with potential V(Q), the 
equation of state varies with the scale-factor, depending on the initial conditions 
on Q, and the form of V(Q). In certain cases, especially for large mass fields 
(e.g. harmonic potential with large mass field, m » Hq), the evolution is 
oscillatory. In most cases, however, the mass of the quintessence field is smaller 
than or comparable to the Hubble parameter (m Hq), and the evolution 
of w(a) is monotonia Observations are consistent with models in which the 
equation of state evolution is monotonic and slowly varying |l3) . 




Figure 1: In the left panel, we see the evolution of the equation of state for a 
set of potentials, all with flQ = 0.6 and w = —1/3 today. The evolution of the 
ratio of the quintessence energy density to the critical energy density is shown 
in the right panel. In both panels, the upper, solid curve represents a constant 
w = —0.55 model, while the lower solid curve represents a constant w — —0.66 
model. The w = —0.55 model has been chosen to best-fit approximate the 
exponential potential, while the w = —0.66 model has been chosen to best-fit 
approximate the quadratic and quartic potentials. Notice the similarity in the 
energy density evolutions. 

In Figure [j] we consider equation of state histories in a group of models with 
different potentials and initial conditions such that w(a) increases monotonically. 
The examples plotted are quadratic, quartic, and exponential potentials ]l9[ . 
We compare the evolution of the cosmological energy density in quintessence, 
Pq(cl), in these models with that from models with constant equation of state. 
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The parameters and initial conditions have been chosen so that each case pro- 
duces the same present-day values for the equation of state and total energy 
density. The left panel shows the evolution of the equation of state, while the 
right panel shows the evolution of the energy density. In both panels, the upper, 
solid curve represents a constant w = —0.55 model, while the lower solid curve 
represents a constant w — —0.66 model. The w = —0.55 model has been chosen 
to best-fit approximate the exponential potential, while the w = —0.66 model 
has been chosen to best-fit approximate the quadratic and quartic potentials. 
For each of the potentials, Pq{o) is a monotonically increasing function of the 
scale-factor. Although the time-history of the equation of state is quite different 
between the constant w and evolving potential cases, we can see from the figure 
that the evolution of the energy densities is closely comparable. 

Quintessence affects the CMB anisotropy chiefly through effects which de- 
pend on pq{cl), that change the expansion history of the universe Thus, 
these potentials predict nearly identical CMB anisotropy to the best-fit con- 
stant w models with w = w e ft. The effective equation of state is empirically 
obtained as the pq weighted average value of w(a) ]T^ |: 

fdaw(a) p Q {a) 
w eff = — ri T~\ — ■ ( 12 ) 




Figure 2: Here we plot the ratio of the power spectrum in the models from 
Figure [j] to the power spectrum in the corresponding best-fit constant w model. 
The fractional cosmic variance with respect to the best-fit model is also shown 
(outer thin lines). The ratio for each model falls well within this variance en- 
velope at most of the multipole moments, thus, the predicted anisotropy is 
observationally indistinguishable. 
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In Figure g we plot the ratio of the multipole moments of the CMB power 
spectrum in the exponential, quartic, and quadratic potential models to the 
power spectrum of the corresponding model with constant equation of state 
w = W e ff. Also plotted is the fractional cosmic variance. (Cosmic variance is 
the intrinsic theoretical uncertainty for any model prediction based on adiabatic 
gaussian perturbations.) The ratio in each of the cases falls within the cosmic 
variance uncertainty for almost all multipolcs. Hence, the power spectrum in 
each case is observationally indistinguishable from the corresponding constant 
w model. 

Thus, for a large class of models with monotonically changing w, the evolu- 
tion of quintessence and its fluctuations are described, to within cosmic variance, 
by a constant effective equation of state. In this paper we restrict ourselves to 
these models, and additionally require the sound speed in the quintessence fluc- 
tuations, or the group velocity of the fluctuations, c 2 to be ~ 1 at sub-horizon 
scales. We do not deal with oscillatory equations of state or models in which 
the kinetic energy is non-canonical and c 2 < 1 at smaller wavelengths, such as 
k-essence models Q. 

4 Solving fluctuation equations for constant w 

To study the evolution of quintessence fluctuations, we numerically evolve the 
fluctuation and Einstein equations in a QCDM model with Q,q = 0.6, h = 0.65, 
and Vlsh 2 = 0.02. In Figure || we plot the quintessence and matter fluctuation 
energy density obtained in this model for a mode with wavelength larger than 
the horizon today (k = 10 Mpc ). The lower three curves are the fluctuation 
evolutions at three different equations of state, w = —1/3, w = —2/3 and 
w = —0.9. We see in the figure that for all the equations of state, the fluctuation 
amplitude first oscillates and decreases, reaches a minimum, and then ultimately 
starts to increase. The decrease of the amplitude is sustained for a longer time, 
and the subsequent increase is sharper, as w becomes more negative, i.e., closer 
to -1. The upper three curves, which are all almost on top of each other, are the 
corresponding matter fluctuation evolutions. The change in the quintessence 
fluctuation evolution as a function of w does not impact the matter fluctuation 
evolution at all. 

To understand the nature of this evolution, we look for analytical solutions to 
the quintessence fluctuation equation. We collect together the other constituents 
of the universe into an adiabatic background fluid denoted by the label l B\ with 
equation of state nig and energy density ps- We can then solve the fluctuation 
equation analytically in two limits: (a) the energy density in the quintessence 
field is negligible compared to that in the other constituents of the universe 
(pq << pb)i arid (b) the energy density in quintessence dominates the energy 
density in the rest of the constituents (pq » ps)- In these limits we can solve 
for 8tj) deep inside the radiation- and matter- dominated epochs where is 
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Figure 3: We compare the evolution of quintessence and matter fluctuation 
energy density for a long wavelength k = 10 Mpc mode in three different 
models with equations of state w = —1/3, w = —2/3, and w = —0.9, and 
with FIq = 0.6, h — 0.65, and f^/i 2 = 0.02. The three lower curves are the 
quintessence fluctuation evolutions at the different equations of state, while the 
three upper curves, all very close to each other, are the corresponding matter 
fluctuation evolutions. Notice that the energy density in quintessence fluctua- 
tions changes with equation of state, but remains much smaller than the energy 
density in matter fluctuations. 



a constant, and deep inside the quintessence-dominated epoch respectively. In 
both cases, the fluctuation equation (Eqn. [ll]) then simplifies to the form of a 
forced harmonic oscillator with constant coefficients and a force term dictated 
by the coupling of quintessence to the perturbed metric: 

d 2 5ip dSip ,, 2 2n X _ 1 dh dip 

d^af + "d(l^a) +{k + m W - - 2 d(lQga) d(lQga) . (13) 

Here k is the co- moving wave number, 

y=l(i-w B ) (14) 

is the damping coefficient of the oscillator equation, and 



m= -y/(l-w)(2 + w + w B ) (15) 

is the Compton mass of the field 5ip, all expressed in units of the co-moving 
Hubble parameter (^-). 
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There are two qualitatively different types of solutions to this fluctuation 
equation. The homogeneous solutions belong to a two-parameter family speci- 
fied by the initial conditions dtpmit and 5ij)' init , and are unaffected by the fluctu- 
ations in the background cosmological fluid. The inhomogeneous solutions, on 
the other hand, arise as the response of the quintessence field to the fluctuations 
of the background. The evolution of quintessence fluctuations as a function of 
scale-factor is determined by the combination of these solutions. The solutions 
are different in the radiation-, matter- and Q-dominated epochs, and must be 
matched by continuity at the boundary between radiation and matter domina- 
tion, and between matter and quintessence domination. 

Solving the inhomogeneous equation requires knowledge of the time evolu- 
tions of the quintessence field and the metric perturbations. The former evolu- 
tion can be obtained from the equation of motion for the field (Eqn. |Io| ) : 



_ / 3i? 2 fJ Q a 2 ( 1 \ I (i+«0 



'(«) ' (16) 



d(loga) V 8ttG 

while the latter dependence can be obtained from the Einstein equations fT^] : 

d 2 h (l-3w B ) dh ^^p + Sp Sp B , , 
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where is the fractional energy density in the background fluctuations. The 

behavior of at scales larger than the Jean's length, the largest scale at 
which the col/apse of a fluctuation through gravitational instability can be coun- 
teracted by the propagation of mechanical disturbances in the baryon-photon 
fluid, in the synchronous gauge, is well known pTJ: 



SpB _ f^PB_\ ( a 
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(18) 



Here (^fjf^j IS the fractional background energy density and a B , is the scale- 
factor when the fluctuation mode under consideration exits the horizon during 
inflation (Hi). The density power spectrum at horizon crossing is taken to be a 
scale invariant Harrison-Zeldovich spectrum with a COBE normalized amplitude 
D B , 

'^-) = D B k^,n = l, (19) 

PB 'Hi 

with p = 4 deep in the radiation dominated epoch, and p = 2 deep in the matter 
dominated one |21fl . 

We use Eqn. ]l^ and the solution of Eqn. |l^ to solve the fluctuation equation 
(Eqn. [K^) for the evolution of metric perturbations, and consequently for the 
inhomogeneous solution in the radiation and matter dominated epochs. At 
wavelengths much larger than the Jean's length, we find that: 

6i>! = - Cl D B a H 2 a% , (20) 



where p = p + 3(wb — w), and 



C/ 6 (r>f,) (v 0* ) w5) m 2 + 2/h/ + 4m 2 • (21) 



Notice that the inhomogeneous solution depends on two separate epochs, an r , 
and clq = 1, the latter entering the equation through the dependence on the 
evolution of the quintessence field. 

The magnitude of the inhomogeneous solution is proportional to the am- 
plitude of the Harrison- Zeldovich spectrum and is thus determined by the 
COBE normalization. Since /i > 0, the solution increases with increasing scale- 
factor for all equations of state. This behavior corresponds to the gravitational 
amplification of large wavelength quintessence fluctuations due to CDM po- 
tential wells. Furthermore, the inhomogeneous solution at a given scale-factor 
is smaller for values of w closer to -1, since the a 2 scaling and coefficient cj 
are both smaller for more negative values of w and the a^ 2 term is almost in- 
dependent of w. The reduced amplitude reflects the smaller coupling to the 
background in the source term of the fluctuation equation (Eq. |^). 

The homogeneous equation has the form of a damped harmonic oscillator 
with constant coefficients. Thus the solution at all wavelengths in each epoch 
is simply: 

SipH = c H a~%9(a,k,m,is), (22) 

where ch is the amplitude of the solution and where 9 is an oscillatory function 
of order unity. Since v > in all epochs, the oscillation envelope decreases 
as a power law of the scale-factor. The amplitude ch must be determined 
by the initial conditions on the quintessence fluctuations at the initial hyper- 
surface far outside the horizon, deep in the radiation dominated epoch (we chose 
0>init ~ 10~ 8 in our simulations and analysis). 

The inhomogeneous solution scales as a positive power of a, and is hence 
negligible (< 10 -20 ) at the initial hyper-surface. The initial fluctuations in 
quintessence are thus entirely due to the homogeneous solution (Sipinit = 5ipH,init) 
Hence, a change in initial conditions affects the homogeneous solution only. To 
determine ch, we consider the initial conditions predicted by inflation. Infla- 
tion creates a nearly scale-invariant primordial spectrum of adiabatic density 
perturbations in all light fields. Since the quintessence fields of interest in this 
paper are also light fields, the entropy perturbation for the entire fluid, just 
after inflation, vanishes: 

TSS = 6p- 2-6p = 0. (23) 
P 

This condition gives one equation between the initial fluctuations S^mu and 
Sip' init . A second constraint is obtained from the observation that long wave- 
length fluctuation modes are frozen outside the horizon, and thus we set: 

Winii = 0. (24) 
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We solve the constraint equations for the amplitude of the homogeneous 
solution: 

6+3™ _ 2 

The declining power law scaling (a - ^) of the homogeneous solution is indepen- 
dent of w. Thus, the equation of state dependence of the homogeneous solution 
comes only from its amplitude ch- Consequently, the value of the homogeneous 
solution at a given scale- factor is larger for w closer to -1. 

Having obtained the approximate solutions of the fluctuation equation as a 



function of scale-factor and equation of state (Eqns. |20| and 22), it is now pos- 
sible to understand the numerically obtained long wavelength evolution shown 
in Figure ||. Firstly, note that the scale of both the solutions is determined by 
the amplitude of the matter fluctuation at horizon re-entry, and consequently 
the COBE normalization. Secondly, we can see from the equations that that 
the homogeneous solution decreases as a~% , while the inhomogeneous solution 
increases as a~ . Thus the amplitude of the fluctuations decreases until the inho- 
mogeneous solution becomes comparable to the homogeneous solution, and then 
it starts to increase. The scale factor at which the solutions become comparable, 



f-c H a Hl \7-6w 

aT = I^t^tJ ' (26) 



increases from ax ~ 10 -5 at w = —1/3 to a*r ~ 1.7 x 10~ 3 at w — —2/3 to 
ot ~ 4 x 10~ 2 at w = —0.9. Thus, as can be seen in the figure, the homogeneous 
solution is comparable to the inhomogeneous one for the w = —1/3 model at 
last scattering (a~7x 10~ 4 ), while it dominates the inhomogeneous solution 
at both the more negative equations of state, w — —2/3 and w = —0.9. 

The magnitude of the homogeneous solution at a given value of the scale 
factor increases as w approaches -1. By contrast, the magnitude of the inhomo- 
geneous solution decreases. Additionally, since 

Sp Q ex y/ oc 1+3m , (27) 
a 2 

the energy density in the homogeneous solution at a given scale-factor further 
increases with decreasing w for all w < —1/3 | fl6| . Consequently, in models 
with w closer to -1, such as w — —0.9, the energy density of the quintessence 
fluctuations is initially larger and decreases more gradually. Thus, the evolution 
remains dominated by the homogeneous solution until a later time. 

The power law decline of the homogeneous solution is independent of wave- 
length. On the other hand, the amplitude of the inhomogeneous solution for 
small wavelength modes is suppressed compared to amplitude for large wave- 
length modes. For wavelengths smaller than the co-moving free streaming scale 
for the quintessence fluid, L f s , the fluctuations free-stream from over-dense to 
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under-dense regions. Thus, modes smaller than L/ s experience oscillations and 
the damping of the power law growth of the inhomogeneous solution due to the 
competing effects of gravitational amplification and pressure support from free 
streaming. 

In this section we obtained approximate analytic solutions to the fluctuation 
equation at long wavelengths. We used these solutions to explain the evolution 
of quintessence fluctuations for different equations of state. We found that for 
w closer to -1, the homogeneous solutions dominate the inhomogeneous ones 
until later in the evolution of the universe. In the next section we study the 
sensitivity of the CMB anisotropy to initial conditions. We show that this longer 
lasting domination at values of the equation of state closer to -1 determines the 
extent to which the initial conditions must be changed from the case of perfectly 
smooth initial conditions to affect the CMB anisotropy. 

5 Sensitivity to Initial Conditions 

We have seen in the last section that initial conditions affect only the homo- 
geneous solutions of the fluctuation equation. For models with w closer to -1 
such as w = —2/3 and w = —0.9, the homogeneous solution is larger and 
dominates the inhomogeneous solution longer. In particular, the homogeneous 
solution dominates at last scattering, and a change in initial conditions can 
propagate forward in time to a change in the total fluctuation energy density, 
and consequently, to a change in the temperature anisotropy. Since the power 
law decline of the homogeneous solution is independent of wavelength, and the 
amplitude of the inhomogeneous solution is suppressed at smaller wavelengths, 
any conclusions on sensitivity to initial conditions drawn at larger wavelengths 
will continue to hold at smaller ones. 

At long wavelengths, an expression for the effect of the fluctuations in p rn 
and pq on the metric perturbation can be obtained in a very simple form from 
the Einstein equations fl~7|| : 

-\?L h > „4*G<?5p m {2+^). (28) 
2 a dp m 

The fluctuations in quintessence produce an effect which depends upon the ratio 
5pQ/8p m . This ratio must become comparable to unity at last scattering for 
there to be any distinguishable effect on the CMB anisotropy. As can be seen 
in Figure |3| in the case of adiabatic initial conditions, quintessence fluctuations 
are sub-dominant to matter fluctuations by many orders of magnitude for all w. 
While the ratio SpQ/5p m will increase if one amplifies the initial conditions, we 
will see that it is still too small at last scattering for most w in order to have 
a distinguishable effect on the metric perturbation and consequently the CMB 
anisotropy. 
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Figure 4: The figure in the left panel shows the evolution of the ratio of energy 
in Q fluctuations to that in matter fluctuations {j^) 2 at fc = lO^Afpc -1 for 
both smooth initial conditions (inhomogeneous solution) and adiabatic initial 
conditions at w = —0.9. The figure in the right panel shows the corresponding 
CMB power spectra as a function of multipole moment. Plotted below the 
power spectrum is the percentage residual of the power spectrum for adiabatic 
initial conditions from smooth ones, compared to the fractional cosmic variance 
(100 x plotted as a black line). The anisotropy change in going from 

smooth to adiabatic initial conditions is well below the variance 



To study the effect of changing initial conditions on the CMB anisotropy 
we start with the simplest possible initial conditions, smooth initial conditions, 
where the values of the fluctuation amplitudes Sip and Sip' are set to zero on 
the initial hyper-surface. Smooth initial conditions have the unique property 
that the quintessence fluctuation evolution is determined solely by the inhomo- 
geneous solution of the fluctuation equation. To test sensitivity we compare to 
the case of adiabatic initial fluctuations in Q. This corresponds to mixing the 
homogeneous solution into the inhomogeneous one. 

In the left panel of Figure 0, we compare the evolution of the ratio 8pQ /8p m 
at large wavelength (k — 10 Mpc~ x ) for smooth and adiabatic initial condi- 
tions, at w = —0.9. We see that the evolution of the ratio for the smooth case 
tracks the power law rise of the inhomogeneous solution to its present-day value. 
The magnitude of the ratio at last scattering (a ~ 10 -3 ) is much larger in the 
adiabatic case than in the smooth case, corresponding to the dominance of the 
homogeneous solution over the inhomogeneous one. Still, 5pQ/Sp m remains far 
below unity in all epochs for the adiabatic case. 

We plot the effect on the CMB anisotropy at w = —0.9 in the right panel 
of the figure, for both the smooth and adiabatic initial conditions. Below the 
spectra is plotted the absolute value of the residual, or the percentage difference 
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Figure 5: The figure compares the evolution of the ratio (y^jj) 2 for both 
adiabatic (F = 1) and artificially amplified (F = 10 4 and F = 10 5 ) initial 
conditions for the w = —0.9 model from Figure ||. The ratio is plotted for 
wave number k = 10~ 4 Mpc~ 1 . We also plot a solid horizontal line to indicate a 
ratio of magnitude unity. Notice that the amplification prolongs the domination 
of the homogeneous solution, and the resultant closeness of the energy density 
ratio to unity. 



of the power spectrum in the model with adiabatic initial conditions compared 
to the model with smooth initial conditions. We also plot the fractional cosmic 
variance. Residuals smaller than the variance cannot be observationally mea- 
sured. As can be seen from the figure, the anisotropy for the adiabatic case at 
w = —0.9 is not observationally distinguishable from the smooth case within 
cosmic variance. While the addition of the homogeneous solution to the inho- 
mogeneous one does increase the ratio 5pQ/5p m by many orders of magnitude, 
the increase is not enough, even at w = —0.9, to alter the CMB power spectrum. 

If there is to be a distinguishable imprint of a change in initial conditions 
on the CMB anisotropy, the energy density in quintessence fluctuations must 
increase drastically so that SpQ / 8p m is of order unity at last scattering. Let us 
estimate how large the initial amplitude of 5pQ / pq must be in order to have 
a distinguishable effect by artificially multiplying adiabatic initial conditions 
by a factor F, and then comparing the result to the results for smooth initial 
conditions. We want to show that F must be quite large in order to have any 
detectable effect on the CMB anisotropy. In Figure || we display the evolution 
of SpQ/Sp m in the w = —0.9 model, for initial conditions that are adiabatic 
(F = 1), and for initial conditions F = 10 4 and F = 10 5 times the adiabatic 
initial conditions. The ratio is plotted at a wavelength longer than the horizon 
today, k = 10~ 4 Mpc -1 . In the case of amplified initial conditions, the homo- 
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geneous solution is larger, and thus it takes until very recent epochs for the 
inhomogeneous solution to become comparable to the homogeneous solution. 
The amplification of the initial conditions by F = 10 4 makes SpQ/Sp m larger 
than unity initially, and the dominance of the homogeneous solution keeps it 
close to, but smaller than unity at last scattering. An amplification by F = 10 5 
makes the ratio to be of order unity at last scattering, which will leave an imprint 
on the CMB anisotropy. 

To understand these effects from the perspective of the value of , we need 
to obtain the value of the matter fluctuation at the initial hyper-surface (a ~ 
10~ 8 ). Since COBE normalization sets the amplitude of the matter fluctuation 

on re-entry to be ^ ~ 1(T 5 , we find from Eqn. pi that (^st) ~ lCT 16 

P m V ? m J init 

in the synchronous gauge. The ratio of ( ^2. ) to ( ^£i~ ) j s se ^ either by 

\ p Q J init \ pm J init 

imposing smooth or adiabatic initial conditions. In the former case we have 
8pq/pq — and so the ratio is zero. In the latter case, the ratio can be 
obtained by combining Eqns. and |24| with the scale-factor dependence of p rn 
and pq. The absolute value of this ratio ranges from ~ 10" 1 at w = to 10 22 
at w = —0.9. In other words, for w closer to -1, ^2 j s initially quite large. Yet 
there is no observable difference in anisotropy, as can be seen from Figure. |], 
between the cases of smooth and adiabatic initial conditions. The amplitude 
of the homogeneous fluctuations swiftly declines and both the above ratio, and 
consequently SpQ/Sp rn are much smaller than one by last scattering. Thus 
there is no observable change in the CMB anisotropy. It is only when the initial 
conditions are amplified by F — 10 5 (so that (j^J ~ 10 11 at the initial hyper- 
surface) that the steep decline cannot offset the initially large value by the epoch 
of last scattering, and there is any observable effect on the CMB. Of course, this 
large value of F is physically unrealistic, many orders of magnitude greater than 
what is expected from inflation, for example. Also, for such extreme values of F, 
the linear approximation used in CMB analysis is invalid. This exercise shows 
clearly that we can ignore the initial conditions on the quintessence fluctuations 
for all reasonable models. 

In Figure ||, we plot the numerically computed power spectra at w — —2/3 
and w = —0.9 for adiabatic and amplified initial conditions. The amplification 
is by factors of F — 10 2 , 10 4 and 10 5 times the adiabatic initial conditions. 
Below the spectra, we plot the residuals with respect to the adiabatic model. 

We see that at w — —2/3, the power spectrum in the amplified models is 
identical to that in the adiabatic model, even for F = 10 5 . For w = —0.9, 
an amplification by a factor of 10 2 leads to no distinguishable changes in the 
anisotropy. On the other hand, an amplification by F = 10 4 weakly suppresses 
the Doppler peak and creates changes in the anisotropy at some multipoles. The 
ratio of energy in quintessence fluctuations to energy in matter fluctuations for 
F = 10 is of the order ~ 10~ 2 -10 -1 at last scattering, as can be seen in Figure 
^|, and the residual anisotropy is barely smaller than the cosmic variance. Thus, 
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w=-2/3 w=-0.9 




multipole moment / 

Figure 6: The figure depicts the CMB power spectrum as a function of multi- 
pole moment for two of the models of Figure || with w = —2/3 and w = —0.9. 
The power spectra are plotted for a series of cases with artificially amplified 
initial conditions, and for the corresponding model with adiabatic initial con- 
ditions. Also shown in the lower panel is the absolute value of the percentage 
residual of the amplified cases from the adiabatic case, as compared to the 
fractional cosmic variance (black line). At w = —2/3, an amplification of the 
adiabatic initial conditions even by F = 10 5 , is not enough to make an observ- 
able change in the CMB power spectrum. On the other hand, in the w = —0.9 
case, the power spectrum for the same value of F is markedly different. 

the effects on the CMB are not large enough to be observationally distinguished. 
Amplification of the initial conditions in the w = —0.9 case by a factor of 10 5 
raises the amplitude of the homogeneous solutions sufficiently that j^- ^ 1 at 
last scattering. For equations of state even closer to -1, smaller amplification 
factors are required to make a measurable difference in the CMB anisotropy. 
However, the amplification is still large and unphysical. For example, even at 
w = —0.999, an amplification by F = 10 3 is required to create an observable 
effect. 

6 Conclusions 

We studied in this paper a large class of quintessence models with light fields 
and sound speed c 2 s ~ 1 at small wavelengths, which have the property that 
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they can be well approximated by constant equation of state, w. The evolution 
of the fluctuations in these models was obtained by numerical integration and 
explained by approximate analytic solutions to the fluctuation equation at large 
wavelengths. Our central result is that the CMB anisotropy in such models is 
insensitive to initial conditions on the quintessence fluctuations for smooth and 
adiabatic initial conditions . For w = —0.9, the CMB anisotropy is insensitive 
in the large range of initial conditions ( ^l) < 10 11 (F = 10 5 ) for ( ~ 

V PQ ) init \ Pm J 

10 -5 at horizon re-entry. Secondly, the sensitivity increases as w approaches 
-1. Atw= -0.999, the range reduces to (^] < 10 9 (F = 10 3 ). However, 

physically reasonable models such as those based on inflation and ekpyrosis do 
not produce such large values of and the ratio of energy in quintessence 
fluctuations to that in matter fluctuations is much smaller than unity. Hence, we 
do not anticipate that the CMB anisotropy will be sensitive to initial conditions 
in realistic cases. The same analytical arguments made in this paper carry 
over to the more general quintessence models in which w is more strongly time 
dependent or c 2 s ^ 1 at small wavelengths. However, the precise numerical lower 
bound on the initial conditions required to imprint a distinguishable effect on 
the CMB anisotropy has to be worked out on a case by case basis. 

This work was supported by DOE grant DE-FG02-95ER40893 (RD), NSF 
grant PHY-0099543 (RRC), and DOE grant DE-FG02-91ER40671 (PJS). 
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